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ABSTRACT 

Most of the proposed associations between magnetars and supernova remnant 
suffer from age problems. Usually, supernova remnants ages are determined from 
an approximation of the Sedov- Taylor phase relation between radius and age, for 
a fixed energy of the explosion ~ 10^^ erg. Those ages do not generally agree 
with the characteristic ages of the (proposed) associated magnetars. We show 
quantitatively that, by taking into account the energy injected on the supernova 
remnant by magnetar spin-down, a faster expansion results, improving matches 
between characteristic ages and supernova remnants ages. However, the mag- 
netar velocities inferred from observations would inviabilize some associations. 
Since characteristic ages may not be good age estimators, their influence on the 
likelihood of the association may not be as important. 

In this work we present simple numerical simulations of supernova remnants 
expansion with internal magnetars, and apply it to the observed objects. A 
short initial spin period, thought to be important for the very generation of the 
magnetic fleld, is also relevant for the modifled expansion of the remnant. We 
next analyze all proposed associations case-by-case, addressing the likelyhood of 
each one, according to this perspective. We consider a larger explosion energy 
and reasses the characteristic age issue, and conclude that about 50% of the 
associations can be true ones, provided SGRs and AXPs are magnetars. 



Subject headings: supernovae remnants - stars: neutron 
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1. 



Introduction 



Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs) are two classes of 
objects being candidates to magnetars, neutron stars with magnetic fields above the quantum 
critical threshold = 4.41 x 10^^ G. In the model proposed by Duncan & Thompson (1992), 
which we will follow through this work, dynamo action in a fast rotating proto-neutron star 
amphfies the "seed" magnetic field (10^^ — 10^^ G) by one to three orders of magnitude. 
The simplest estimate of the magnetic field intensity is derived from measurements of their 
rotational periods (P) and period derivatives (P), through the well known expression 



where B is the magnetic field intensity, Rpsr is the pulsar radius (~ 10® cm), x is the angle 
between the rotational and magnetic field axis, / is the moment of inertia (~ W^^ g cm^), 
and c is the speed of light. Their ages can also be estimated through the characteristic age 



but both calculations assume continuous spin-down driven by magnetic dipole braking, an 
assumption which has been questioned by several authors (Marsden et al. 2001; Harding, 
Contopoulos & Kazanas 1999; van Paradijs, Taam & van den Heuvel 1995; Chaterjee, Hern- 
quist & Narayan 2000; Kouvehotou et al. 1999), in different contexts, regarding AXPs and 
SGRs. The derived figures are shown in Table 1. As most of the objects are younger than 
10^ years, it seems reasonable to look for the Supernova Remnant (SNR) that was originated 
at the same supernova explosion. In fact, almost all SGRs and AXPs were tentatively asso- 
ciated to some SNR. Though, the likelihood of those associations is put in doubt because the 
age estimated for those SNRs is generally larger than the characteristic age of the associated 
magnetar. 

This situation leads to the "age" and "velocity" problems of SGRs and AXPs. Both 
problems are related to the relatively low characteristic ages found for these objects. The 
SNRs usually associated to them seem to be systematically older, thus the age problem. 
Some SGRs and AXPs have been found to lie beyond, at or near the border of their proposed 
SNRs, and, given their low ages, the implied velocities are much higher than those found for 
ordinary pulsars (see Hurley 1999 for a discussion about these objects). As seen on Tables 
1 and 2, there is a large uncertainty in those estimates, mainly due to uncertainties in the 
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distances. This implies uncertainties in the radii on which age estimates are based through 
models of the expansion. 

Several authors have addressed magnetar-SNR associations. Some of them have dis- 
carded either the magnetar hypothesis or the association, based on those problems. This 
became known as the "nature versus nurture" debate, where one side proposes that the 
characteristics of SGRs and AXPs are derived from their special nature (see Thompson & 
Duncan 1996), while the other side finds more plausible that the environment is the ac- 
tual responsible, nurturing the otherwise ordinary neutron stars with a fossil accretion disk 
(Marsdcn ct al. 2001; Chaterjee, Hernquist & Narayan 2000; van Paradijs, Taam & van den 
Heuvel 1995). 

Associations between pulsars and SNRs are generally examined according to the fol- 
lowing criteria: 1) positional coincidence in sky; 2) distance estimates; 3) age estimates; 4) 
evidences for interaction between the neutron star and the SNR. Additional criteria may 
be derived from these, like the estimate of the projected velocity on sky of the compact 
object. Gaensler et al. (2001) did a comprehensive study of positional coincidences for 
magnetar-SNR associations. Distance can be determined by neutral hydrogen column den- 
sity measurements, evidence of association with another object (e. g., HII region or molecular 
cloud) whose distance is already known, or through surface brightness-diameter (E — D) re- 
lations. This last method is not considered reliable enough, because it ignores the effect 
of density and structure of the local interstellar medium, resulting in a large dispersion of 
values from the best fit, as can be seen in Case & Battacharya (1998). Some SNRs on the 
proposed associations show plerions, usually taken as positive signal of the presence of a 
neutron star, but most do not present that feature. Finally, measurements of HI column 
density are subject to systematic errors from models, and usually are not better than 50%, 
implying on similar uncertainties in distance determinations. Radio pulsars can have their 
distance estimated by use of dispersion measure, however SGRs and AXPs have not yet been 
detected in radio waves. There are reasons to believe they cannot emmit radio signals, either 
because they are magnetars, where the huge magnetic field induces photon splitting rather 
than pair creation (Baring & Harding 1998, but see Camilo et al. 2000), or in the case they 
are accretors, because accretion quenches radio emission (Chaterjee, Hernquist & Narayan 
2000; van Paradijs, Taam & van den Heuvel 1995). 

SNR age estimates are usually derived from the knowledge of the radius by means of 
an approximate version of the expansive blast-wave evolution equations. We shall see in 
next Section that those approximations which ignore the previous phase (s) can introduce an 
apreciable error in cases near phase transition, and such approximations imply a reduction of 
the number of variables to two by linking the explosion energy, the local interstellar medium 
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density and the ejected mass. The present local density of the interstellar medium can be 
estimated within 50% uncertainty, though this may not be representative of the structure of 
the interstellar medium prior to the supernova explosion, for example, an explosion which 
occurred inside a HII region, (low density medium) and have run into a nearby molecular 
cloud (high density medium) would be misleading, since the observer would be able to 
estimate the molecular cloud density only. Even inside the HII region it would be likely that 
density is not constant. Situations like those can be modelled (see a good review on this in 
the work of Truelove & McKee 1999), but the point is that the unknown structure of the 
pre-supernova environment introduces one more possible systematic error on age estimates. 
The ejected mass can be estimated from the observation of how much Fe was produced, 
coupled to models of supernova yields, when the SNR is young. Explosion energy is usually 
considered to be within a factor of 2 of the canonical value of 10^^ erg, due to previous both 
observational and theoretical knowledge. The review of Hamuy (2003) considers the energy 
range 0.5 — 5.5 x lO^"*^ erg for classical SNII, from a sample with a dozen well-observed events. 
However, these estimates also contribute with their own uncertainties, because the type of 
explosion event giving rise to a SGR-AXP is completely unknown. As the radius picks the 
same relative uncertainty as the distance, the combined uncertainty on the age estimate will 
not be smaller than a factor of 2, and probably much larger than that. We must recall that 
it is not uncommon to be unsure about the phase a given SNR is in, and that near the 
transition time between phases approximations do not strictly hold. More accurate relations 
can be used, reducing systematic errors, though usually most authors feel the rather small 
gain in precision is not worth the trouble. 

Blast-wave velocity estimates could also be used to determinate age and radius, through 
relations similar to the ones for radius and age. Knowledge of both radius and velocity would 
determine the phase of the expansion, removing a source of uncertainties. Unfortunately, 
such estimates must come from X-rays measurements difficult to be made, and up to now 
they could not be obtained for any of the proposed associations. 

Finally there are evidences of interaction between the neutron star and the SNR in some 
of the proposed associations, which show plerions (or filled-centre morphology). Additionally, 
one of the associations was reported as presenting a jet-like feature (Kouvehotou et al. 1998). 

The overall reliability of the described estimates is not good, although several papers 
have addressed the association issue, and drawn conclusions about AXP and SGR nature. In 
this work, we will show that at least an important factor for those estimates is still missing, 
namely the injection of energy by the internal magnetar, and we analyze the proposed 
associations bearing that factor in mind. The next Section is dedicated to SNR expansion 
considerations. Section 3 introduces the effects of the injected energy on the SNR, while 
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Section 4 deals with an application of those considerations to the observed sample. In 
Section 5 we discuss our results. 



2. Supernova Remnant Expansion 

Supernova explosions eject several solar masses to the interstellar medium, which form 
the SNR. In the further evolution, mass from the swept interstellar medium will be added, 
at a rate AnR'^Rp, where R is the SNR radius (blast- wave front), R is the SNR expansion 
velocity, and p represents the local interstellar medium density, which we will take as a 
constant for simplicity. A treatment including density as a power-law function of the distance 
from the supernova can be seen on Truelove & McKee (1999), which is also a good review 
reference about the non-radiative expansion phases of SNRs. The initial expansion velocity 
is set by the total kinetic energy of the explosion, E, and the ejected mass, M^j. The 
interstellar medium will not affect much the expansion until the swept mass becomes ~ M^j, 
or, equivalently, the SNR mass 



M = Mej + —R^p ~ 2Mej . (3) 

3 

Another approximation that holds true is that very little energy (compared to the huge 
kinetic energy of the SNR) is lost, so energy can be considered constant. In this way, 
with mass and energy set to constants, and negligible external pressure, velocity is also 
(approximately) constant. This is the free expansion phase. An approximate linear relation 
between radius and age {t) follows 



A more accurate expression can be found in the work of Truelove & McKee (1999): 



R ~ 0.46 pc eH^M;^'\ i^l + 0.011 El{^M;^'HTn\'^) , (5) 

where n = pm'^p~^ is the number density of the interstellar medium, mn is the hydrogen 
mass and p — lAis the mean molecular weight of the interstellar gas. Energy is scaled in 
units of 10^^ erg {E^i = E/IQ^^ erg), mass in units of 10 Mq (Mio = M/10 Mq), number 
density in units of 1 cm~^ (rii = n/1 cm~^), and age in units of 100 years {t2 — i/10^ yr). 
The velocity at which the blast-wave front expands is approximately 



-6- 



R 4500 km s'' eH'm-,"' (l + 0.011 Eli'M-.'/STn^') ■ (6) 

As M increases, so does the ram pressure of the interstellar medium, slowing down 
the SNR expansion. Growing mass and pressure must then be taken into account. The 
SNR enters the Sedov-Taylor phase, named after the works of Sedov (1959) and Taylor 
(1950) about pressure- driven explosions. To describe the SNR expansion one should solve 
the equation 

^(^Mr) = AnR-'p, (7) 



dt\4: 



the internal pressure being 



where 7 is the adiabatic index (5/3 for an ideal gas, 4/3 for a relativistic one) and U is 
the internal energy of the gas inside the internal cavity formed by the remnant, whose mass 
concentrates in a thin shell (its thickness will be neglected). As the total energy is roughly 
constant, one can use 

U = E - —MR^ . (9) 

Remembering that eq. 3 shows M depends on Mej and R, eq. (7) turns to have no ana- 
lytic solution. Nonetheless, as at that phase M^j <C M, one could use those approximations 
to obtain the well-known analytic solution 

R ^ (C^) t^i' , (10) 

where C, = 2.02 in the exact solution, and we found 1.77 (both for 7 = 5/3; if 7 = 4/3, 
our solution would be 1.06). By construction, this solution disregarded the previous phase. 
Truelove & McKee (1999) have calculated a corrected (though still approximate) expression, 
which is 



~ 12.5 pcMio^V^^^ ^5(^^17' 0.051 . (11) 



.1/2, ,-5/6^1/3 



2/5 
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where age was given in units of 10^ years (^4 = yr). The correspondent expression for 

velocity is 



it ~ 490 km s-^ eH^M'^''^ 



El{^M~^'^n\'^U - 0.051 



-3/5 



(12) 



Although we have been considering the constancy of the total energy, in fact the SNR 
has been slowly radiating away. As the SNR expands, its temperature decreases, because it 
depends on the blast-wave velocity, following the well-known relation for strong shocks. 



16 Kbu 



— -^f2.02-V^-«/^ 
100 kBn\ n) 



(13) 



where fc^ is Boltzmann's constant, and R was taken from the derivative of eq. (10). 

Eventually T will reach 10^ K, where the dominating cooling process changes from 
thermal bremsstrahlung to line emission, which is more efficient to radiate away energy. 
That way, the adiabatic approximation ceases to be accurate, and the SNR enters into the 
snowplow phase. As a rough approximation to estimate when the Sedov- Taylor phase ends, 
it is usual to estimate how much time it would take to reduce the thermal energy of the SNR 
to zero, considering just radiative losses. Starting from the energy loss per particle. 



-nA , 



(14) 



where A = 1.6 x 10~^^C^~^/^ erg cm^ s~^ is a simple cooling function appropriate for ionized 
gas at temperatures 10^ K < T < 10'^ K, = 1 is a mctallicity factor parametrized for solar 
abundances, one can integrate that equation and equate the result to the dynamical time 
R/R ~ 5t/2, obtaining 



3/14 -4/7 



tsT-SP ^ nxWyvE^^in, 



(15) 



If we choose to use eq. (12) into eq. (2), a slightly more accurate expression results 



tsT-SP ^ 19 X 10^ yr eII'\ + 510 yr E^l'^ M%\-^^'^ . (16) 

It has been suggested that SNRs are very difficult to be seen after 20 kyr (Braun, Goss 
& Lyne 1989). The similarity of the values of "fading time" and the transition from the 
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Sedov-Taylor to snowplow phase may be taken as an indicative that they are related by a 
factor, of order unity. Thus, wc will adopt eq. (16) as the rough limit of visibility of one 
SNR. Since our intention is to address associations between SNRs and magnetars, we will 
not explore further the expansion of SNRs into the snowplow phase, because no SNR would 
be detectable in that phase or afterwards. 

3. Energy injection by a Magnetcir 

The energy loss of a pulsar is usually taken as arising from a rotating magnetic dipole 
approximation, 

L.sr- ^ 

That energy is intercepted by the SNR, which will contain, in an internal cavity, most of 
either relativistic particles or electromagnetic waves emmited by the central object. Either 
way, this cavity would thus be filled by a relativistic gas, pushing the SNR from the inside. 

Considering that the magnetic field, moment of inertia and x ^i-re constants, they can be 
absorbed together with other factors into a new constant K — ^ ^gjj'" ^ ■ After integrating 
eq. (17) the period evolution can be expressed as 

P = ^JsTT^Kt + Pg , (18) 
with Pq being the initial period of the pulsar. Two additional constants will be defined 

^0 = = O.Gdays ( Y -^V , (19) 

is the initial time-scale for deceleration, and 

U^K,[^y - 3.85xl0«e.gs-' (^)^(^)", (20) 

is the initial rate of energy loss. We notice that Lqt^ oc B~'^ and LqTq oc Pq"^. With these 
relations and eq. (18) we can rewrite eq. (17) as 
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Lp,, = Lo(^l + -j . (21) 

Total injected energy is just the integral of eq. (21) from the initial instant to 1 s 
< To to the present time, 



-1 ^ -1 



2.0 X 10^^ erg 



1 ms 



B 

lO^Gj \lmsj Vlyr 



1.6xl0-3( 1 1^1(^1 +1 



(22) 



Eq. (19) and (20) arc already scaled to the typical magnetar parameters. We remark 
that unless magnetars are born with very short periods (~ 1 ms), their magnetic fields would 
not be expected to grow enough to cross the critical quantum boundary Be, according to 
the model of magnetar formation of Duncan & Thompson (1992). While pulsars can be 
born with such short periods, they are not required to do so. The difference of two orders of 
magnitude in the magnetic field strenght between a typical pulsar (~ 10^^ G) and a typical 
magnetar (~ 10^^ G) means 4 orders of magnitude in both tq and Lq values. This implies 
in turn a dramatically different influence on SNR expansion: a magnetar will inject most of 
its rotational energy into the internal cavity of the SNR within a day, while a pulsar would 
take tens or hundreds of years, depending on its initial period and magnetic field. Also, a 
magnetar will inject typically 10^ times more energy than a pulsar, and, more remarkably, 
that energy is a factor of 10-20 bigger than the kinetic energy of an ordinary SNR. It is 
analogous to the suggestion of Ostriker & Gunn (1971) about that transfer of energy being 
the very cause of a supernova event. 

Gravitational radiation losses are usually larger than rotating magnetic dipole ones 
for pulsars with very short periods, hence any estimate of the initial period based on the 
electromagnetic dipole torque will result in a figure corresponding roughly to the period of 
transition between gravitational and magnetic dipole dominances. In fact, a simple estimate 
can be obtained equating eq. (17) to the following expression for gravitational wave-carried 
energy loss, taken from Shapiro & Teukolsky (1983): 

^arav - g ^5 I p i ' '^^'^^ 

where G is the familiar gravitational constant and e is the oblateness of the neutron star. 
The result is 
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Ptran = 16a/^ — = 0.34 HIS ( ^VV^) , (24) 

V 5 RpsrcB sin X yW'^G j VlO-V ' ^ ^ 

which shows that for a magnetar the large magnetic dipole losses would not be much affected 
by gravitational losses. In other words, the initial period will effectively be the same period 
that has allowed the magnetic field intensity to increase above Be, and almost all rotational 
energy will be injected on the SNR. 

However, there is room to consider competition for the energy loss dominance between 
magnetodipole radiation and gravitational radiation by r-modes. Doing a very simple anal- 
ysis, from the work of Owen ct al. (1998). wc would expect r-modes gravitational waves to 
be more efficient than magnetodipole radiation to extract rotational energy from the rapidly 
rotating magnetar for P < 4 ms. Though we feel more detailed studies of both the possible 
damping of r-modes by the neutron star crust, and the enhancement of magnetic fields in 
newly-born neutron stars are in need before we can address properly this issue. 

Earlier works have addressed pulsar energy injection on a SNR (Luz & Berry 1999; van 
der Swaluw et al. 2001), but they have not explored the case of a magnetar source. It can 
be easily seen that the scenario built by van der Swaluw et al. (2001) cannot hold when 
magnetars are considered, since they have assumed the energy injected by the pulsar to be 
smaller than the kinetic energy of the SNR. 

The introduction of magnetar-injected energy changes the equations describing the ex- 
pansion of the SNR in the following ways: eq. (8) remains the same, but 7 = 4/3 could be 
used to represent the dominance of the injected energy over the initial kinetic energy of the 
SNR; eq. (9) picks up a new term, becoming 



U = E - -MB? + ^ ° , . (25) 

However, in this simple formulation of the problem, only the blast-wave front is de- 
scribed, and we have not described the full problem of two gases. The energy injection is 
so quick that we consider all the SNR is instantaneously reacting to it. That is not strictly 
true, but since the internal shock (the inner cavity boundary) would be near the external 
shock (the SNR boundary) in less than 100 years, that is much before the transition to 
Sedov-Taylor phase starts, this simplification will not affect much Sedov- Taylor and poste- 
rior phases. We will only consider 7 = 5/3 through this work. While the works of Luz & 
Berry (1999) and van der Swaluw et al. (2001) studied the injection of energy by a pulsar 
into a SNR, we adopted simplifications that are not well suited to be combined with their 
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methods and results. However, the general results obtained in this work are not expected to 
be affected by these approximations. 

We have performed numerical simulations to find solutions for the set of equations 
describing the position and velocity of the blast- wave front, with and without including 
the energy injected by an internal magnetar, the former case to test the simulation engine. 
Our results for the case with energy injection are similar to the ones in the case without 
energy injection, though with explosion energy set to the initial kinetic energy plus energy 
injected (from the internal magnetar case). Those results are nearly identical because of the 
extremely short time needed for the injected energy to raise above the initial kinetic energy. 
Our numerical results differ only from the expressions shown on eqs. (11) to (16) in one 
point, namely our values for R and R in the free expansion phase differ from Truelove & 
McKee (1999) as if the energy were reduced by a factor of 1.14. This factor comes from the 
approximation done on Section 2, specifically the difference in ^ values on eq. (10). 

This factor must be multiplied to the energy if our values are to be compared to Truelove 
& McKee (1999). Through this work, however, we will adopt the values we have found, 
without corrections. The discrepancy is not large, and certainly does not affect much our 
considerations, given other simplifications done in this work. 

The numerical results can be appreciated in Figs. 1 and 2. In Fig. 1 we show the 
evolution of R{t) for both cases (with and without energy injection). In Fig. 2 we show the 
evolution of R{t) for both cases. 

4. Analysis of Proposed Associations 

The proposed associations between would-be magnetars and SNRs are shown in Tables 
1 and 2. We shall analyze them by looking first into the situation without considering the 
injection of energy by a magnetar. This situation is shown in Fig. 3, where we show the 
range of radius and ages associated to each SNR considered, the characteristic ages of the 
magnetars, and the radius evolution curves for two different cases (thick solid lines): a "low 
density/low mass" evolution scenario, in which M^j = 8 Mq and n = 0.01 cm~^, and a 
"high density/high mass" evolution scenario, with M^j = 30 Mq and n = 10 cm~^. For 
both extreme cases, explosion energy is held fixed to 10^^ erg. The position of a given 
SNR should be between these two extremes, unless quite different explosion energies are 
considered. Certain SNRs have no reliable age estimate (Marsden et al. 2001); for those, 
we adopted an arbitrary range 0.2-30 kyr, which is the possible range for galactic SNRs. 
For several associations, SNR and magnetar ages are not compatible, and some magnetar 
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ages are not compatible with the SNR expansion in enviroments with ordinary values for n. 
That has been used to justify the model of AXPs and SGRs being born in regions of higher 
density than radio pulsars (Marsden et al. 2001), or to dismiss the association altogether. 

Varying the ejected mass value will displace those curves diagonally (in free expansion 
phase only) , but will not affect meaningfully the position of the curves in the Sedov- Taylor 
phase. In any case, quite unreasonable values of both Mgj and n should be invoked to 
maintain some associations as valid. The other way to displace the curves (in both phases) 
is by changing explosion energy. 

In Fig. 3 wc also show the curves including the energy injection by a magnetar with 
i? = 5 X 10^"^ G and Pq = 1 ms. The displacement of curves helps to attribute lower values 
to the density than before to all associations, and it makes possible to "save" some otherwise 
untenable associations. Thus, the injection of energy could be behind the age discrepancy, as 
SNRs truly associated with magnetars would have expanded faster than expected. In other 
words, the SNR which has been born with an internal magnetar will seem older. The actual 
relation between the true age {tt) and the apparent one {to) can be obtained from eq. (11), 
for the Sedov- Taylor phase, considered for both the conventional energy value {Eg) and the 
one including energy injected {Et), as 




Typical figures would be in the range ~ 0.2 — 0.3 1^. 

Likewise, other reference timescalcs will be shifted. The transition from free expansion 
to Sedov- Taylor phases will occur sooner, because of the higher initial velocity. The transition 
to snowplow phase will occur at 36 kyr after the supernova explosion (see eq. 16), using 
Et — 20Ea. This is just meant to show that an SNR with energy injection could be visible for 
more time than an ordinary SNR, while at the same time appearing to be younger. The time 
a neutron star will take to catch up its SNR (tcross) before the transition to the snowplow 
phase (while it is still visible) can be roughly estimated from eq. (11) and the neutron star 
(constant) velocity v — 10^V2 km s~^, resulting in 

tcross = 6.5 X lO^yr v-^^^EK^n-'^^ (27) 

where wc ignored the small negative term on eq. (11) for the sake of simplicity. From this 
last equation it can be seen that only a fast magnetar can catch up its SNR shell. For 
example, if V2 = 10, as suggested for some associations, and E^i = 20, then tcross — 38 kyr. 
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Of course, the proximity of the neutron star to the SNR shell can induce a "reenergization" 
of the shell, extending the visibility of the SNR to later ages (see ShuU, Fesen & Saken 1989). 
We do not address this possibility here. 

Writing the distance of the neutron star to the center of the SNR as r{t) — l3{t)R{t), 
the quantity /3 becomes an observable parameter (projected in the sky). Inverting eq. (27) 
and inserting /3, we can find an expression for the minimal neutron star velocity needed to 
reach a relative displacement /? at age t^, 

^; = 1.2 X 10^ km s"^ pt^^^^EH^'^^^ . (28) 

If the neutron star velocity points transversally to the line of sight, then eq. (28) gives the 
actual velocity of the neutron star. 

4.1. Analysis of magnetcir candidates with known P 

We shall analyse the proposed associations taking into account the energy injected 
by the internal magnetar, addressing the plausibility of the association, and the nature of 
the compact object. Similar studies have been published by several authors (Mereghetti 
1999; Marsden et al. 2001; Gaensler et al. 2001), based mainly on position and velocity 
considerations. Those works have disregarded magnetar characteristic ages, in favour of 
SNR estimated ages. However, their conclusions are discrepant. We shall in turn disregard 
the ages of the SNRs, because if they were born with an internal magnetar, their ages are 
overestimated according to eq.(26). 

The novel feature of our approach is just the fact that, if the AXPs and SGRs are mag- 
netars, following the model of Duncan & Thompson (1992) for magnetar formation, which 
is based on a dynamo action to make possible the growth of the magnetic field beyond ordi- 
nary pulsars range, they will inject energy enough on their SNRs to affect their expansions, 
making them reach a larger size in less time than it was considered by previous authors. Al- 
though our simulation engine is probably too crude to provide a reliable estimate of the time 
after which the injection of energy becomes less efficient, even more sophisticated simulations 
have not been tested to know this number. We will take as enough to say that injection 
times below one month are securely within the bounds to provide optimal coupling between 
injected energy and the SNR kinetic energy. As the main effect of the energy injected is to 
add to the initial kinetic energy of the SNR, any pulsar which can inject most of its rotational 
energy in less than one month will provide an SNR evolution as if the kinetic energy were the 
initial plus the injected one. The injection time for the smallest B of the AXP/SGR sample 
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(AXP 2259+586) is lower than 5 days, considering Pq ~ 1 ms. Even initial periods as large 
as 3 ms would provide injection times close or within to the one month figure, which make 
us certain of the eff'ectiveness of the coupling. Pulsars, on the other hand, have injection 
times longer than decades, typically. 

For simplicity, and having no way to estimate the actual initial periods of AXPs and 
SGRs, we assume all them as being born with Pq = 1 ms, which in turn determines that 
the injected energy (~ 2 x 10^^ erg at t » tq) will be essentially the same, regardless of 
the actual value of B (from cq. 22), as long as Tq < 1 year. Because of this fact, we will 
only consider two extreme situations to analyze the associations: the "low density/low mass" 
evolution scenario, where M^j = 8 M0 and n = 0.01 cm^^, and the "high density/high mass" 
evolution scenario, with M^j = 30 Mq and n = 10 cm~^. The actual situation for each case 
should be bracketed between these values. The result can be seen in Fig. 4, in which the 
range of radius and characteristic ages values for each association is explicitely shown. Our 
assesment of each association is as follows: 

SGR 1806-20/G10.0-0.3: The probability of ahgnment by chance is ~ 0.5% (Marsden 
et al. 2001). GlO.0-0.3 was considered as not being a SNR (Gaensler et al. 2001; Chakrabarty 
et al. 2001), although Green (1998) lists it in his catalogue. A cluster of stars is close to 
the line of sight to this association (Fuchs et al. 1999), so either the SGR or the SNR may 
be physically related to it. In our model, considering the characteristic age of SGR 1806- 
20 as the true age of the association, it can be seen that the entire range of radius values 
lie between the two extreme scenarios (Fig. 4), with an indication of mid to low-density 
interstellar medium. The magnetar transversal velocity implied is high, 4000-6500 km s^^, 
if /3 = 0.5 (Kulkarni et al. 1994). However, if /3 ~ (Ghakrabarty et al. 2001), the velocity 
cannot be inferred. GlO.0-0.3 would be entering Sedov-Taylor phase. Dropping altogether 
the characteristic age as a good age estimator, the association could be as old as 15 kyr, 
with V ~ 500 km s~^ {P = 0.5). This association can be considered as true if SGR 1806-20 
is a magnetar, ~ 0, and GlO.0-0.3 is confirmed as a SNR. 

SGR 1900-I-14/G42.8-I-0.6: The probability of random alignment is ~ 4% (Gaensler 
et al. 2001). A pulsar was recently discovered near that position (Lorimer & Xilouris 2000), 
that could be related to the SNR, although its characteristic age is 38 kyr. Our model 
would allow for the association if the true distance were on the lower 1/3 of the range and 
the true age were on the upper 1/2 of the quoted range, but the extremely high velocity 
implied (> 8000 km s~^) precludes that. Disregarding characteristic age, the range would 
be extended up to 6-40 kyr, depending on true distance and interstellar medium density, 
though the magnetar velocity would still be high (800-2000 km s~^). Wang, Li & Zhao (2002) 
suggested that SGR 1900-1-14 was born in the 4 B. G. supernova, and the discrepancy between 
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characteristic age and the proposed age of 2 kyr is attributted to dynamical evolution with 
braking index ~2, but they did not offer a good reason to explain how the SNR disappeared 
from view in just 2 kyr. We conclude that this association is not convincing (if SGR 1900+14 
is a magnetar), unless a mechanism for a high velocity of the magnetar is adopted. 

AXP 1048-5937/G287.8-0.5: The probability of chance ahgnment is ~ 16% (Mars- 
den et al. 2001; Gaensler et al. 2001). Data about G287.8-0.5 was considered unreliable by 
Gaensler et al. (2001). Nevertheless, if the SNR is confirmed as such, our model indicates 
that it is entering the Sedov- Taylor phase, on a high density interstellar medium. Low ages 
arc preferred. Again, the very high velocities implied (4500-8500 km s~^) argue against 
the association. In this case, even to disregard the characteristic age does not improve the 
plausibihty of the association. We conclude that this association is unhkely. 

AXP J1709-4009/G346.6-0.2: The random alignment probability is ~ 10% (Mars- 
den et al. 2001) to ~ 30% (Gaensler et al. 2001). Our model could allow the association if 
the interstellar medium has high density, and the true distance is on the upper 25% values of 
the range. Once more, the high velocity required for /3 = 1.7 (3000-2000 km s~^) is in excess 
of the known pulsar population. Disregarding the characteristic age would only help very 
slightly to avoid the "velocity problem" . We therefore consider this association as unhkely. 

AXP J1709-4009/G346.5-0.1: Gaensler et al. (2001) suggest this association instead 
of the previous one. While the newly-identified SNR G346. 5-0.1 awaits to be confirmed, 
we can analyze the association in the same fashion we have been doing. The probability 
of random ahgnment is ~ 10% (Gaensler et al. 2001). Our model allow the association if 

the true distance is on the upper 80% values of the range. Once more, the high velocity 
required for f3 = 1.2 (1700-4800 km s^^) exceeds the typical pulsar velocity. Disregarding 
the characteristic age would allow the age range to go up to ~ 50 kyr (~ 800 km s^^), at the 
transition to snowplow phase. We consider this association as more likely, requiring a mid to 
high-density medium. We notice that our solutions with injection are not better than the old 
situation (without energy injection), so the association can be true even if AXP J1709-4009 
is not a magnetar. 

AXP 1841-045/G27.4+0.0: The chance alignment probability is ~ 0.01% (Gaensler 
et al. 2001). In our model, this association would require an excedingly low density inter- 
stellar medium, or negligible energy injection. Equivalently the association would require 
Po > 6 ms, or age< 600 yr (but with v > 1600 km s~^). Even then, the characteristic age 
would have to be disregarded. Contrary to most AXPs and SGRs, the characteristic age 
of AXP 1841-045 is higher than the SNR estimated age. The association is possible only if 
AXP 1841-045 is not a magnetar with short Pq. 
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AXP 2259+586/CTB 109: The probability of random alignment is ~ 0.05% (Gaensler 
et al. 2001). The characteristic age of AXP 2259+586 puts its associated SNR in the snow- 
plow phase, where it would be unlikely to be detected. Given the low probability of chance 
alignment, this age discrepancy (this is the other AXP which has characteristic age higher 
than SNR age) argues again against the characteristic age as a good estimate of the true 
age. Removing this parameter allows for our model a wide range of possible age values, 
from 1 to 30 kyr. The requirement of reasonable magnetar velocity would limit ages < 6 
kyr. Moreover, this range is already allowed by the expansion without energy injection. The 
association is probably true for both models. 



4.2. Analysis of magneteir candidates with unknown P 

For these objects, we have no information about their characteristic ages or magnetic 
field strenght. Therefore the task here is to verify if the allowed range of ages of our model 
is compatible with reasonable magnetar velocities. 

SGR 0526-66/N 49: The random alignment probability is ~ 0.7% (Gaensler et al. 
2001). The allowed range of ages is 0.5-3.5 kyr, which implies extremely high velocities, 

V > 1800 km s~^. If this association is true, then SGR 0526-66 would have to be another 
magnetar born with a long Pq, as suggested for AXP 1841-045. 

AXP J1845-0258/G29. 6+0.1: The chance alignment probability is 0.2% (Gaensler 
et al. 2001). The allowed range of ages is 1-13 kyr, which imphes on 3500 > v > 300 km 
s~^, for the largest distance of the range, and 0.4-1.5 kyr (3600 > v > 1000 km s~^) for the 
opposite extreme of the distance range. This association can be true with or without energy 
injection. If AXP J1845-0258 is a magnetar, then G29. 9+0.1 prefers the higher end of values 
for distance and age. 

SGR 1627-41/G337.0-0.1: The probability of random alignment is ~ 5% (Gaensler 
et al. 2001). The allowed age range of our model to this association (150-300 yr) imply on 

V > 15000 km s~^, ruling out the magnetar hypothesis, unless Pq > 6 ms, as discussed for 
AXP 1841-045. Even so, the velocity problem still holds. Thus, we find the association 
unlikely. 

SGR 1801-23/G6.4-0.1: SGR 1801-23 is just a candidate SGR, and its position is 
not well-determined. The distance to G6. 4-0.1 is uncertain, too. Within such set of data, 
our model allows for wide ranges of ages. In the high extreme of the distance range, ages 
can be 2-35 kyr (1000 > f > 70 km s~^). In the low end of the distance range, ages allowed 
are 0.45-2 kyr (1500 > v > 350 km s~^). Thus, the association is probable, and this holds 
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even if SGR 1801-23 is not a magnetar at all. Though, this result is heavily dependent on a 
better position determination. 

5. Discussion 

We have analyzed in this work the proposed SGR-AXP/SNR associations. Previous 
anahsis of the associations arrived at different results. According to Marsden et al. (2001), 
all associations can be considered likely. Gaensler et al. (2001) contend that only AXP 
J1845-0258/G29.6+0.1, AXP 1841-045/G27.4+0.0 and AXP 2259+586/CTB 109 could be 
valid. Ankay et al. (2001) considered SGR 1806-20/G10.0-0.3 and SGR 0526-66/N 49 as 
plausible, in addition to the 3 already mentioned by Gaensler et al. (2001). 

It is expected that if magnetars exist and are born in supernova explosions, they can 
inject enough energy to enhance the expansion of their associated SNRs. If AXPs and SGRs 
are indeed magnetars, their associated SNRs should appear older than the ages derived from 
standard expansion models. We analyzed all the proposed associations, and have come to 
the following results: 

1) If the characteristic age of the neutron star is regarded as the true age of the as- 
sociation, then SGRs and AXPs may not be magnetars at all (or the model of Duncan & 
Thompson (1992) is not correct regarding magnetar origin), for only one case (SGR 1806- 
20/G10.0-0.3) has shown good agreement within the model, and even the true nature of 
GlO.0-0.3 was put in doubt; 

2) If characteristic ages are in turn ignored (sec Harding, Contopoulos & Kazanas 1999; 
Gaensler ct al. 2001; Marsden et al. 2001; Mcrcghetti et al. 2002, among others), then 
two SGRs (out of 5) and three AXPs (out of 5 with proposed associations) seem to be 
associated within our model, so about 50% of general agreement. It should be noted that 
two associations which agreed within our model (AXP J1845-0258/G29. 6-1-0.1 and AXP 
2259-1-586/CTB 109) were already believed to be true by previous works (Gaensler et al. 
2001; Marsden et al. 2001; Ankay et al. 2001), meaning that uncertainties in distance are 
large enough to allow for both possible scenarios (standard and with energy injection); 

3) AXP 1841-045/G27.4+0.0 and SGR 0526-66/N 49 can only be considered as true 
associations if the magnetars were born with Pq > 6 ms, because in that way the energy 
injected would be insufficient to directly affect SNR expansion. 

These results from our model are tied to the dynamical evolution of magnetars. While 
we have assumed for simplicity the standard magnetic dipole braking with braking index 
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equal to 3, several proposals have been done that argue for different braking models: fossil 
or fallback accretion disks (Marsden et al. 2001; Chaterjee, Hernquist & Narayan 2000; van 
Paradijs, Taam & van den Heuvel 1995), episodes of relativistic wind emission (Harding, 
Contopoulos & Kazanas 1999), a different constant braking index (Wang, Li & Zhao 2002), 
magnetic field decrease (Duncan & Thompson 1992; Colpi, Geppert & Page 2000). As the 
measurements of braking indices for 5 young pulsars revealed that all but one have significant 
departures from the canonical value, it should not come as a surprise that magnetars do not 
follow the standard model of spindown. The main effects of alternative models are to change 
the estimates of magnetic field strength and spindown age. Nonetheless, the infiuence in our 
model would be very small, since Tq can be increased up to one month without appreciable 
modifications on the model, and rotational energy is not dependent on spindown models. 
That is why we feel justified to leave item (1) above aside of the discussion. Spindown models 
that allow for ages significantly smaller or larger than the conventional characteristic age are 
required to explain certain associations, if they are true ones, either within our model, or 
considering the standard scenario. 

Supposing that our results from item (2) above represent the actual situation, we could 
ask where are the SNR that are associated to the other SGRs and AXPs. The answer can 
be that they were not detected yet, since they lie generally at regions of coexistence of HII 
regions, other SNRs, and other types of objects (variable stars, young stars, molecular clouds, 
etc.). Including AXP 0142+614 and AXP J0720-312 in the sample, we would have slightly 
less than 50% of the true associations identified. On the other hand, few of the latter can be 
considered as firm, due to several other factors. For example, the angular size of SNRs are 
often quoted without uncertainties, which certainly is an understatement, given the dificulty 
to recognize a SNR and to assess its shape and size. 

It is worthy to remark that a similar analysis, ignoring the effects of the injection 
of energy, and desconsidering conventional characteristic ages, can be done to the same 
associations, and will result in different age and velocity ranges than those obtained here. 
Indeed, even associations considered unlikely by us may turn out to be likely. This is because 
the standard energy SNRs would take more time to reach the observed sizes, allowing for 
higher ages for the associations than in our model, and thus producing lower velocities than 
those we found. Nevertheless, that scenario requires non-standard spindown and no injection 
of energy by the neutron star, be it a magnetar or not. 

It is interesting to check that the most important factor to decide on the plausibility 
of associations is /?. As /5 increases, associations are considered increasingly unlikely. How- 
ever, Gvaramadze (2002) points to the possibihty that a SNR expanding in a region with 
anisotropic interstellar medium densities will be distorted and/or expands faster in one or 
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more directions. This means that the geometrical center of the SNR can be displaced from 
the actual explosion site. Although it is difficult to take this effect into account quanti- 
tatively, one should be cautious to do not dismiss it entirely. It is possible that one or 
more of the proposed associations arc affected by this effect, which can increase or decrease 
f3 randomly. That effect will be more accentuated for SNRs that have been expanding in 
low-density regions ("bubbles") amidst high-density "walls", and it increases with age. The 
first consideration is the case for the most massive stars, which do not live long enough to 
move far away from the sites where they had been born, while their stellars winds create 
a low-density cavity. The second consideration means we must expect /3 estimates to be a 
little bit more scattered as age (and radius) increases. 

Besides studying associations including AXPs and SGRs, we also examined the associ- 
ation PSR J1846-0258/G29.7-0.3 (Gotthelf et al. 2000; Mereghetti et al. 2002), which also 
suffers from the age problem, since it has the smallest characteristic age known among pul- 
sars (723 years), and the SNR age was estimated as at least 1800 years. The magnetic field 
of this pulsar is ~ 5 x 10^'^ G, considering magnetic braking spindown, slightly above the 
quantum critical field, and by so marginally qualifying as a magnetar. It was not detected 
at radio wavelenghts, only in X-rays. The association is considered very likely, as the pulsar 
is located at the geometrical center of the SNR, coincident with a radio/X-ray nebula, prob- 
ably powered by the pulsar. Proceeding as in the other cases investigated in the previous 
Sections, we see from Fig. 4 that the range of radius quoted on the hterature is nearly 
coincident with the range allowed by our two extreme cases, considering the characteristic 
age. This way, there are no preferences for high or low densities or ejected mass. As before, 
we will ignore the characteristic age, to found that the age range allowed is 250-700 years, 
in case of high M^j or n, and 650-6000 years, in case of low Mej or n. Intermediate values 
of the parameters M^j and n would provide intermediate ranges of ages. The placement of 
the neutron star at the geometrical center of the SNR imply in low velocities or alignment 
between the velocity vector and the line of sight. It is important to notice that Mereghetti 
et al. (2002) find a braking index ~ 1, 9 and age ~ 1700 years for this pulsar, in which 
case the magnetic field can be under the quantum critical one. Nonetheless, this association 
represents one more evidence against the consideration of the characteristic age (without 
braking index information) as a good age estimate. 

We have left for a future work (Allen & Horvath, in preparation) the study of an 
alternative origin scenario for magnetars, the collapse of a white dwarf star, induced by 
accretion or merging from a binary companion. Simulations for this scenario (Fryer et. al. 
1999) reveal that ~ 0.1 Mq can be ejected, with an explosion energy of 10^° erg, implying 
on higher initial velocities. 
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Finally we would like to point that magnetars lose very little rotational energy through 
gravitational waves, when compared to typical pulsars, unless r-modes can play an important 
role, which is not clear to us. Statistical studies on pulsar gravitational wave dctcctability, 
like the one performed by Regimbau & de Frcitas Pacheco (2001), shall not be affected by 
this consideration, at least for the current detectors, like VIRGO, because of the scarcity of 
magnetars. 

Wc arc grateful to G. Medina- Tanco for discussions about SNR blast-wave physics. This 
work was supported by FAPESP Agency (Sao Paulo State, Brazil) and CNPq (Brazil). 
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Table 1. Data of SGRs and AXPs 



magnetar P (s) P (IQ-^^-) ^ ^^^^^^ ^ (^^qm ^ ^j^p^-) 



SGRs 












SGR 1806-20 


7.48 ' 


83 ^ 


1.4 * 


8.0 ' 


17 ^ 14 ' 


SGR 1900+14 


5.16 ^ 


50 - 140 9 


0.58- 1.6 1 


5.1 - 8.6 « 


7 ^ 5 


SGR 0526-66 


8.00 








55 ^ 


SGR 1627-41 


6.4 ? ^ 








11^ 


SGR 1801-23 












AXPs 












AXP 1048-5937 


6.45 ™ 


15 - 40 « 


2.6 - 6.8 « 


5.1 - 8.6 


10.6 > 2.8 


AXP J1709-4009 


11 


19 *, 22.5 ^ 


9.2 *, 7.7 ^ 


4.6 *, 5.0 ^' 


10 > 8 « 


AXP 1841-045 


11.8 


47 41 <^ 


3.9 4.5 


7.6 ^, 7.0 




AXP 2259+586 


6.98 


0.74 0.49 * 


150 230 * 


0.73 0.59 * 


6.2 " 


AXP J1845-0258 


6.97 ^ 








8.5^ 


other 












PSR J1846-0258 


0.324 ^ 


7.1 ^'^ 


0.72 ^'^ 


0.5 ^'^ 





References. — (a) Duncan & Thompson (1992); (c) Kouvehotou et al. (1994); (d) 
Gotthelf & Vasisht (1997); (h) Vasisht & Gotthelf (1997); (i) Kouvehotou et al (1998); 
(j) Kouvehotou et aL (1999); (k) Hurley et al (1999); (m) Chakrabarty et al (2001); (n) 
van Paradijs, Taam & van den Heuvel (1995); (o) Gaensler (2000); (p) Kaspi (2000); (q) 
Mereghetti (1999); (t) Kaspi, Chakrabarti & Steinberger (1999); (v) Israel et al. (1999); 
(A) Corbel et al. (1999); (B) Woods ct al. (1999); (E) Torn et al. (1998); (X) Gotthelf et 
al. (2000); (Y) Mereghetti et al. (2002). 
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Table 2. Data of associated SNRs 



SNR ta (kyr) d (kpc) R (pc) /3 
associated to SGRs 

G 10.0-0.3 0.2 - 30 ^ 13 - 16 ^ 14 - 19 ^ 0.0 ™, 0.6 

G 42.8+0.6 0.2 - 30 3 - 9 11 - 31 ^ 1.2 1.4 ^'^ 

N 49 5 - 16 55 ^, 50 7.1 - 8 8.5 ■'^ 0.6 ™, 1.0 ^ 

G 337.0-0.1 0.2 - 30 ' 11 ^'"'^ 2.4 1.7 ^ 2.3 

G 6.4-0.1 > 2.4 ^ 1.2 - 3 ^ 3.5 - 4 ^ 7 - 17.5 ^ 21 - 24 ^ 0.1 

associated to AXPs 

G 287.8-0.5 0.2 - 30'' 2.5-2.8'' 9.1-10.2^ 2.2'' 

G 346.6-0.2 0.2 - 30 '' 3 to 5 '\ 11 ' 4.4 to 7.3 '' 1.7 



r 



G 27.4+0.0 2 ^ < 3 « 6 to 7.5 ^'"^ 3.5 - 4.4 '^'^ 4.7 ^ <0.25 0.1 ^ 

G 109.1-1.0 3 to 20 « 4 to 5.6 16 - 24 ^ 18 - 25 0.3 0.2 ^ <0.2 

G 29.6+0.1 0.2 - 30 ^ 8.5 < 20 6.5 to 9.8 ^ 0.1 ^ < 0.25 ^ 

associated to other 

G 29.7-0.3 1.8 - 7 ^ 9 - 21 " 3.9 - 10.7 "'^ 0.0 ^ 



References. — (b) Kulkarni et al. (1994); (d) Gotthelf & Vasisht (1997); (e) Green 
(1998); (f) Matthewson et al. (1983); (g) ShuU (1983); (h) Vasisht & Gotthelf (1997); 
(j) Kouvehotou et al. (1999); (k) Hurley et al. (1999); (m) Chakrabarty et al. (2001); 
(o) Gaensler (2000); (q) Mereghetti (1999); (r) Marsden et al. (2001); (s) Gaensler et al. 
(2001); (w) Sanbonmatsu & Helfand (1992); (x) Gaensler, Gotthelf & Vasisht (1999); (y) 
Rho & Petrc (1997); (z) Vancura ct al. (1992); (A) Corbel ct al. (1999); (D) Corbel et 
al. (1997); (X) Gotthelf et al. (2000); (Y) Mereghetti et al. (2002). 
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age (years) 

Fig. 1. — Radius evolution for a supernova with Mej = 10 Mq, n = 1 cm~'^, and E = 10^^ erg. 
The sohd hne represents the case without energy injection, and the dotted hne represents 
the case with energy injection by a magnetar with B — 5 x 10^^ G and Pq — 1 ms. Squares 
mark the transition to Sedov- Taylor phase, and diamonds mark the transition to snowplow 
phase, where these curves are no longer valid. 
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Fig. 2. — Expansion velocity evolution for a supernova with Mgj — 10 Mq, n — 1 cm~^, and 
E — 10^^ erg. The solid hne represents the case without energy injection, and the dotted 
line represents the case with energy injection by a magnetar with B — 5x 10^^ G and Pq — 1 
ms. 
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Fig. 3. — Comparison of age and radius ranges for the proposed associations to the models 
with and without energy injection. The thick sohd hues represent radius evolution for the two 
standard extreme cases ("low density/low mass" and "high density/high mass" scenarios, 
see text). The dotted thin lines represent the evolution with energy injection for the two 
extreme cases. The thick dashed line is the approximate end of Sedov- Taylor phase for the 
model without energy injection. SGRs (filled squares) and AXPs (filled diamonds) with 
estimated P values are placed according to the median of their characteristic age ranges 
and the median of the radius ranges for the associated SNRs. SGRs (open squares) and 
AXPs (open diamonds) without estimated ages are placed on the median values for their 
associated SNRs. SNRs are small dots placed on the median values of the ranges (shown as 
error bars). For SNRs with unreliable ages, we assumed an arbitrary range 0.2-30 Kyr. The 
new association proposed for AXP 1709-4009 is marked as a filled circle. PSR J1846-0258 is 
marked as a triangle. 
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Fig. 4. — Comparison of radius and ages of the proposed associations to the expansion model 
with energy injection, for objects with estimated P. The sohd hnes represent the evolution 
according to our model. Large markers are placed according to the median of characteristic 
ages and radii ranges. Small markers are placed according to the median of estimated 
SNR ages and radii. Error bars indicate ranges of characteristic ages and SNR radii. The 
dotted line is the approximate end of Sedov- Taylor phase (including energy injection). Filled 
markers represent likely associations, and open markers represent the unlikely ones (see text). 
The markera are: open triangle (SGR 1900+14), filled triangle (SGR 1806-20), open circle 
(AXP 1048-5937), open square (AXP 1709-4009), open diamond (AXP 1841-045), filled 
diamond (AXP 2259+586), filled square (AXP 1709-4009, new association), filled circle 
(PSR J1846-0258). 
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Fig. 5. — Comparison of radius and ages of the proposed associations to the expansion 
model with energy injection, for objects without estimated P. The solid lines represent the 
evolution according to our model. Markers are placed according to the median of estimated 
SNR ages and radii. Error bars indicate ranges of SNR ages and radii. The dotted line 
is the approximate end of Sedov- Taylor phase (including energy injection). Filled markers 
represent likely associations, and open markers represent the unlikely ones (see text). The 
markers are: open triangle (SGR 1627-41), open circle (SGR 0526-66), filled diamond (AXP 
1845-0258), filled square (SGR 1627-41). 



